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Abstract 



We calculate the electromagnetic pion form factor at intermediate space- 
like momentum transfer from the QCD sum rule for the correlation function 
of two pseudoscalar interpolating fields and the electromagnetic current. This 
correlator receives essential contributions from direct (i.e. small-scale) instan- 
tons, which we evaluate under the assumption of an instanton size distribution 
consistent with instanton liquid and lattice simulations. The resulting form 
factor is in good agreement both with the sum rule based on the axial-current 
correlator and with experiment. 
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A central goal of strong-interaction physics is the understanding of hadrons on the basis 
of quantum chromodynamics (QCD). It is very unlikely that this goal can be reached without 
a thorough understanding of the QCD vacuum. At present, however, direct links between 
observed hadron properties and the vacuum structure are still rare and rely mostly on exten- 
sive numerical simulations. The aim of this letter is to study one such link - between direct 
instantons, i.e. small-scale topological vacuum fields, and electromagnetic pion properties - 
analytically in the framework of a QCD sum rule [HJ^. 

Due to the Goldstone nature of the pion, sum rule calculations of pionic properties 
can be based on two in principle (but not in practice!) equally suitable sets of correlation 
functions, corresponding to the use of pseudoscalar or axial vector interpolating fields. The 
pion couples strongly to both of these source currents and thus contributes to the correlators 
in both channels. 

The pseudoscalar channel has some principal advantages for sum rule calculations. The 
accuracy of the standard pole-continuum parametrization for the corresponding spectral 
functions profits from the almost complete dominance of the pion in the low-mass region 
0,^. Furthermore, correlators involving the pseudoscalar interpolators have a simpler tensor 
structure. This simplifies in particular the calculation of three-point functions. 

All existing sum rule applications in the pion channel (with the exception of ref. 0, 
see below), however, are based on axial vector current correlators The use of 

the pseudoscalar interpolating field has been avoided, since it is known to receive essential 
contributions from direct instantons 0]. Like instanton contributions in general, they could 
initially not be estimated reliably, due to insufficient knowledge of the instanton size distribu- 
tion in the vacuum. The attempt of an ab initio description in the dilute gas approximation 
1^, in particular, failed for all but the smallest instanton sizes due to infrared problems with 
large instantons. The same problems were encountered in the attempt to supplement the 
conventional operator product expansion (OPE) in QCD sum rules with direct instanton 
contributions and led to the preference for the axial vector correlators. 
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In the last decade, however, mutually consistent information about the instanton size 
distribution has been gathered from new sources, including phenomenological estimates 0, 



variational studies |10] and numerical simulations of instanton liquid models as well as 



lattice calculations |T^. The bulk features of the density n{p) of instantons with size p in the 
vacuum, which emerged from these studies, can be summarized in a simple parametrization 
i, 

n{p) = n6{p- p) , (1) 

with the average (anti-) instanton density and size fixed at 

n = - fm""^, p = -fm . (2) 

2 ' 3 ^ ' 

This form provides a reasonable approximation to the sharply peaked gaussian distribution 
found in ref. and should be sufficiently accurate for our purpose. 



Armed with the bulk features of the instanton density, we can now try to turn the 
former vice of the pseudoscalar interpolating field - its strong coupling to instantons - into 
a virtue, by using pseudoscalar correlators as a tool for "instanton diagnostics". Linked 
by the corresponding sum rules to observable hadron properties, these correlators provide 
a theoretical laboratory both for the study of instanton effects at short and intermediate 
distances and for tests of the approximations used to calculate them. They also allow an 
independent check of the instanton distribution (|l|). 

With this motivation in mind, we calculate in the present paper the electromagnetic 
pion form factor from a QCD sum rule based on pseudoscalar interpolating fields. The 
instanton contributions are evaluated semiclassically in the zero-mode approximation. Our 
investigation complements similar studies of the nucleon [13] and pion |^ mass sum rules. 



The sum rule techniques for the calculation of hadron form factors at intermediate mo- 
mentum transfers were developed in refs. 0,01 • Their application to the pion form factor, 
based on axial vector interpolators, leads to a rather good agreement with experiment (in 
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the momentum transfer region 0.5 < < 3.0 GeV^) |J3]. Our calculation will thus also 
provide a consistency check between the axial vector and pseudoscalar sum rules. 

To derive the sum rule for the pion form factor, we consider the three-point correlation 
function of two pseudoscalar currents j5{x) = d{x)i'~^^u{x) and the electromagnetic current 
fi[x) = euu{x)'y^u{x) + edd{x)'^^,d{x), 

T,{p^p'-, q) = -jd'xj d'ye^^P'-^-'^-y\0\ T jI{x) fj {y) ^0) |0) , (3) 

with q = p' — p. can be decomposed into two independent Lorentz vector structures, 

T^{p,p'; q) = V^{p\p'\ g2) (p' + + V^{p\p'\ q^) q, . (4) 

The invariant amplitudes Fi 2 have a double dispersion representation of the form 

^i{P,P,q) = — ds ds- -— ^ + ... , (5) 

TT^ JO JO [S — p'^j^S' — p ) 

where the ellipsis represents subtraction polynomials in p^ and p'^, which will vanish under 
the Borel transform to be applied later. 

Since the pseudoscalar current has a nonvanishing matrix element between the vacuum 
and the one-pion state. 
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(0|rf^75M|7r+) = , K = ^ (6) 

if-K = 93 MeV is the pion decay constant and mu,d are the up and down current quark 
masses), the pion contribution to the spectral density is of the form 

(0|j5V(pO)(vr(pOlj?k(p))(vr(p)|j5|0) = 2^2ir2g.F.(g2)(p' + pV , (7) 

where = —q"^, Q-^ is the charge of the pion and F^^IQ"^) is the form factor to be calculated. 
The latter is contained in the amplitude Fi, on which we will concentrate in the following. 
For the continuum contribution we adopt the standard form of refs. PJ^l' which completes 
our parametrization of the spectral density pi: 

pi(s, s', Q2) = 2 /2 F^(g2) 5(5') + + s'- So) Po(s, s', Q^). (8) 



{po{s, s', Q"^) is the free-quark spectral function and we neglected the small pion mass.) The 
continuum threshold sq defines a triangular region in the (s, s') plane and is related to the 
threshold Si of the corresponding pseudoscalar two-point correlator as sq — 1.5 Si (see ref. 
0] for details). 

An alternative continuum ansatz with a quadratic integration region bounded by Si was 
also considered in for the axial-vector three-point function. Similar to ref. we find 



almost no difference between the two parametrizations in the resulting form factor |ll4 |. 

The leading instanton (and anti-instanton) contributions to T ^ are calculated in semi- 
classical approximation, which amounts to evaluating the correlator (|^) in the background 
of the (anti-) instanton fields and to averaging their collective coordinates over the corre- 
sponding vacuum distributions. We will outline only the essential features of this procedure 
and refer to the literature P,P|,P!B| for more details. 

Since the sum rule probes the correlation function mainly at distances x,?/ ~ 0.2 fm, 
which are small compared to the average inter-instanton separation ^ ~ 1 fm, the correlated 
contributions from two or more instantons should be small. We are thus led to treat only 
the effects of the (anti-) instanton closest to x and y explicitly and to include effects of 
interactions with other vacuum fields (including other instantons) at the mean-field level 



To obtain the explicit instanton contribution, we recall that the quark spectrum in the 



background of an (anti-) instanton contains one zero-mode state per flavor |16], 



W) = - f 2^ 2X3/2 (9) 

with r = X — z, where z is the instanton position. (The spin-color matrix U satisfies 
((T + r) f/ = 0.) Up to corrections of order m*p (see below) from higher- lying states, quarks 
in these two zero-modes dominate the instanton contribution to the correlator (|^). Their con- 
tribution is evaluated by inserting the explicit expression (^ into (^) and by approximating 
quarks propagating in the higher-lying (continuum) states as non-interacting. 
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For quarks propagating in zero-mode states the interaction with the other vacuum fields 
leads at the mean- field level to the generation of an effective quark mass m*(p) = — |vr^p^(gg) 



15| , which replaces the current mass in the zero-mode propagator. Taking this effect into 



account, we obtain (after a Wick-rotation to euclidean space-time) the following instanton 
contribution to the correlator (H): 



y d'^x J d'^y j d^ze'P'-''e-''^-ye'^P'-''^-' x 



^6 772*2 (^p) 

+ - iy + zfy^^ 



+ 



X,. 



(10) 

In this expression we summed over instanton and anti-instanton parts and integrated over 
their positions. Since we deal with a gauge-invariant correlator, the averaging over the color 
orientations is trivial, and the average over the remaining collective coordinate, the instanton 
size p, was weighted with the distribution (|I]). 

After taking the standard double Borel transform 0, both in and p'^, of ([10|) (details 
of the rather lengthy calculation will be given in ]14[); we obtain 



(M is the Borel mass) in terms of the dimensionless integral 



(11) 



-^insti^Z , Z ) 



da 
Jo 



ae 



^H{a')H{^') 
-{2ae-A) 



a + e ( ez^ z"^ 

— 17- [a z^ + 



A^{1 - ez^) 
-8{l-ez^ 



X 



a^e 



16 (l-ez2)2 
ae{a + e) 



/i(a')^(7') + -^HiaVih') , (12) 



A ' Az^A\ 

where A = + ae, a' = 7' = and H{x) = Ii{x) — /o(x) is defined in terms of the 
modified Bessel functions Inix). 



The lowest-order perturbative contribution to Fi was evaluated in ref. 00, 



^Note that the corresponding expression in p| contains a misprint. 
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167r2 



V^r\Q\M') = '-^IMQ'^M') , (13) 



where the continuum contributions from the parametrization (H) are subtracted in the inte- 
gral 

As aheady noted in ref. the power corrections to r^^'^'^*'' are small in the range under 
consideration! and will be neglected in this letter. A more complete analysis of the OPE can 
be found in [Q. 

Adding instanton ( pT]) and perturbative (|13|) contributions and equating them to the 
Borel transformed dispersion representation (^) leads to the sum rule 



. (15) 



The factor L^^^^, with L = ln(M^/AQQj-,)/ln(yU^/AQQj-,), accounts for the scaling behavior 
of the three-point correlator due to the anomalous dimension of the pseudoscalar currents, 
and sets their renormalization point to /x = 500 MeV. (We use the value Aqcd = 150 MeV 
for the QCD scale parameter.) 

Let us now fix the remaining parameters in the sum rule. As already mentioned, the 
continuum threshold sq can be related to the corresponding threshold si of the pseudoscalar 
two-point correlator, sq ~ 1.5 Si. The analysis of (see also ref. 0) finds Si ~ 2.0 GeV^ 
and we will thus fix the threshold at Sq = 3.0 GeV^. The relatively large separation of the 
continuum from the lowest resonance in the pseudoscalar channel (about twice of that in 
the vector and axial vector channels) is clearly favorable for the sum rule analysis, since it 
improves the accuracy of the parametrization (H). 

Due to uncertainties in the determination of the light quark masses, the phenomenological 
value of the (scale-dependent) mass parameter K{A), defined in eq. (^, is at present not 



^See also the related discussion for the two-point correlator in S. 
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accurately known. From the current bounds on the up and down quark masses and with 
= 138 MeV we obtain 1 < K{1 GeV) < 2 GeV. In order to allow for a direct comparison 
with the instanton analysis of the two-point correlator |p, which takes K = 0.7 GeV, we 
fix K at the lower limit of the phenomenologically acceptable range, = 1.0 GeV. The 
effective mass m*, finally, is determined as in ||5|Jl3[] from the self-consistency relation P] for 
the quark condensate, (qq) = —2n/m*{p). 

Let us now proceed to the quantitative analysis of the sum rule ([l5|). Figure 1 shows 
the dependence of both the perturbative and the instanton contributions to the pion 
form factor at fixed = 1 GeV^. Two features of the instanton part are immediately 
apparent: First, and as expected, it is the dominant contribution, about a factor of two 
larger than the perturbative part. Secondly, it improves the stability of the sum rule, i.e. it 
reduces the dependence of the form factor, and a "stability plateau" begins to develop 
for M2 > 1.2 GeV^ The same qualitative behavior is found for all values of > 0.5 GeV . 

In order to determine the dependence of the form factor, we follow the procedure of 
ref. and evaluate the sum rule for different values of at a fixed value of the Borel mass, 
= 1.2 GeV^ postponing a more accurate analysis to ref. [jl3 . The resulting form factor is 



shown in fig. 2 and compared with the experimental data of |18| at the space-like momentum 
transfers accessible within our approach!. The agreement is clearly satisfactory and at least 
as good as the one from the axial vector sum rule (dashed line), which we show for 
comparison. (A slightly better agreement between the axial vector sum rule result and the 
data was obtained in ref. by extracting the dependence of at = 1.8 GeV^.) 

The somewhat better fit of the pseudoscalar sum rule result at low (where also the 
data become more reliable) is perhaps not accidental. The result of the axial-vector based 



^The perturbative part dominates the instanton contributions at — > oo and shows the expected 
(Q'^)~^ behavior. The so far neglected power corrections, however, will eventually blow up |14|, 
thereby limiting the applicability of the sum rule at large Q^. 
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sum rule relies exclusively on the OPE, which breaks down at small Q^, since physics of 
longer distance scales begins to determine the behavior of the three-point function. The 
instanton contribution, in contrast, which dominates our result, could optimistically remain 
reliable up to distances not too far below the inter-instanton separation, corresponding 
to ~ 0.1 — 0.2 GeV^. Indications supporting this expectation have been found in an 
analogous study of the nucleon correlation functions [|l^]. The applicabilty of the sum rule 
at small would then be mainly limited by the unphysical singularity in the perturbative 
contribution ([131) . 

To conclude, we view the above results, and in particular the stability of the pseudoscalar 
sum rule and its agreement with phenomenology, as further support for the importance of 
the instanton component in the QCD vacuum, for the bulk features of their distribution as 
given in eq. (^, and for the semiclassical estimate of their effects in hadron correlators at 
intermediate distances. 

H.F. acknowledges support from the Department of Energy under Grant No. DE-FG05~ 
93ER-40762 and M.N. acknowledges support from FAPESP Brazil. 
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FIG. 1. Borel mass dependence of the pion form factor from Eq. (|l5| ) at = 1 GeV^ (solid 
curve). The dashed and dotted curves correspond to the instanton and the perturbative contribu- 
tions, respectively. 

FIG. 2. The electromagnetic pion from factor from Eq. (|l5|) (solid hue) at = 1.2 GeV^. The 
dashed line shows for comparison the result of the axial vector sum rule [^j. The experimental 
data are taken from ref. [18|. 
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